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Abstract 

We investigate the role of the cross decay rates in a system composed by two electromagnetic 
modes interacting with the same reservoir. Two feasible experiments sensitive to the magnitudes 
and phases of these rates are described. We show that if the cross decay rates are appreciable there 
are states less exposed to decoherence and dissipation, and in limit situations a decoherence free 
subspace appears. 
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It has become increasingly important to understand how the environment acts over a 
system we are interested in. The promise of Quantum Computation often has as its main 
obstacle the entanglement of the variables of interest with the reservoir degrees of freedom 
[3| . Also, in the investigation of Foundations of Quantum Mechanics, the environment plays 
a central role |2|. 

A very successful model for the environment is the Caldeira-Legget model, where one har- 
monic oscillator linearly coupled to a bath of oscillators is considered jsj . The investigation 
of bipartite quantum systems interacting with the environment usually displays amazing 
effects jj, sj]. In 0, the Caldeira-Legget model was extended for two oscillators subjected 
to the same bath, leading to a possible communication between the oscillators mediated by 
the reservoir. This extension may be used to study two electromagnetic modes constructed 
in one or two cavities. 

In pi , an experiment dealing with two modes constructed in a single cavity was described. 
The experimental findings were compared to theoretical results in [a], with satisfactory 
agreement. However, due to the large detuning of the modes, the experiment is not sensitive 
to the cross decay rates: the parameters related to the their communication through the 
environment. 

Cross decay rates and their interference effects over quantum systems are well-known for 
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a long time jjj. Some results related to them are population trapping |5|, 
free subspaces (DFS) 0,13]. The extension of the Caldeira-Legget model for two oscillators 
predicts the possibility of DFS if the cross decay rates are large enough [8| . Given the growing 
interest in DFS, the knowledge about this topic is accordingly important. Thus, we propose, 
here, feasible experiments to investigate these rates in the cavity Quantum Electrodynamics 
context, and also identify states that are more resistent against decoherence and dissipation 
when the cross decay rates are not zero. 

The Hamiltonian we use to model two electromagnetic modes subjected to the same 
environment is 

H = H + Hint, (1) 

H = M^ajai + M2 2 a|a 2 + ft^cjfccj.c fc , 
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where ai, a2, aj and a\ are annihilation and creation bosonic operators for the modes of 
interest, Mi and M2, with frequencies Q\ and SI2- The annihilation and creation operators 
Cfc and cj. are used to model the environment by a set of harmonic oscillators linearly coupled 
to Mi and M 2 modes. This Hamiltonian shall be used for two modes in different cavities or 
in the same cavity. Under the usual approximations, Hamiltonian leads to the master 
equation 

j t p s {t) = C Ps (t), (2) 
with the Liouvillian superoperator (an operator which acts on operators) 
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for zero temperature, where 
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and we used the conventional notation for superoperators [13| (the dot sign (•) indicates 
the place to be occupied by ps(t), where the superoperator acts). A discussion about 
the Hamiltonian (0) and a detailed derivation of the master equation (j2J) from may be 
found in j^. The coupling between the environmental modes and the modes of interest 
may occur by complicated processes, e. g., a photon may be scattered from a mode of 
interest to an environmental mode by an atom. Of course, H int does not come from a first 
principles analysis, and, thus, we do not know much about the coupling constants a\k and 
a.2k- However, as we have said, this model presents results that agree with experimental 
ones. The cross decay rates are ki 2 + «A 12 and k 2 \ + «A 2 i, and will assume relevant values 
if the summations involving aifcC^fc an d «2fc«ifc are n °t nu U. This demands that the ot\k 
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and the a 2 k must have some correlation, i. e., the way the modes of interest interact with 
the environment must be microscopically correlated. This may be achieved if the modes of 
interest are close to each other in the scale of the wavelength of the environmental modes 
that effectively interacts with them. These are the ones with frequencies around fli or Q 2 , 
as may be seen in Eq. (J1J). If the environmental modes are electromagnetic modes, the 
modes of interest must be close in the scale of their proper wavelengths. An interesting case 
is two modes with orthogonal polarizations constructed in the same cavity. They occupy 
essentially the same positions in space, but maybe the difference in the polarizations spoils 
the correlation. 

In Fig. ^ we sketch an experiment where circular Rydberg atoms A s and A p , with 
relevant levels e and g, are produced in box B, cross the cavities Ci and C2 and are 
detected in detector D. The energy of level e is higher than the energy of level g by hQ a . 
The cavities are macroscopically identical: Qi = Q 2 = Q an d k%% = k 22 = k. We will assume 
a huge number of environmental modes, with random distribution in the frequencies around 
Q, what leads to An = A 22 ~ and k\ 2 + iA i2 ~ (k 2 i + iA 2 i)*, as may be seen using 
expression (J3J). We may adjust the atom-fields interaction time using the Stark effect: when 
an appropriate voltage is applied across the mirrors of the cavities, the e <-> g transition 
becomes resonant with the modes Mi and M 2 ; when no voltage is applied, the detuning 
A = Q a — Q turns negligible the atom-fields interaction. 

In the analysis below, global phases will be often disregarded. The cavities are initially 
in vacuum state, and A s is sent in state e: 

l^ = 0)> Mll M 2 ,A s = |0,0,e>. 
Atom A s enters cavity Ci, and is put in resonance with mode M\ by a time ti s , yielding 

\*l> (t = tu)) M lt M 2 ,A s = cos ( Gt ^) |0, 0, e) - i sin (Gt ls ) 1 1, 0, g) , 
where we used the RWA (rotating wave approximation) Hamiltonian 

H ls = Ml (ajai + aJ,a 2 + + (aJo"_ + aicr + j , 

with a z = \e) (e| — \g) (g\, a_ = \g) (e\, a + = |e) (g\. During a time t 0s , A s is put far of 
resonance with the modes. Considering |A| 3> G, the Hamiltonian may be written 

H 0s = Ml ^ajax + a 2 a 2 j + ; 
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and the system evolves to 

|V> (t = h s + t 0s )) MlM2jAs = cos (Gt ls ) |0, 0, e) - j e <(n«-«)*o. sin (Gtls ) | 1; , <?) . 

Next, A s enters C2, and is put in resonance with M 2 by t 2s = it / {2G). Considering that 
Mi and M 2 have the same polarizations, the Hamiltonian will be 

H 2s = (ajai + a|a 2 + ^) + ^ ( a 2°"- + a 20+) , 

what leads to 

|V (t = tis + tos + hs)) Ml ,M 2 = cos # l°- !) + e i0 sin# |1,0> , (5) 

where = Gt ls and = (O a — fi) t 0s . The A s state ends up in g, factorized, and needs not 
to be considered anymore. 

During the time T (T 3> t ls +t 0s + t 2s ), no atom interacts with the modes. Using Eq. (J2J) 
to compute the action of the environment in this period, we get, for the state of M\ <-> M 2 
system, 

p s (t = T + t u + t 0s + t 2s ) = {ui (T) |1, 0) + u 2 (T) |0, 1)} {HC} (6) 

+ |0,0) (0,0|{l-| Ul (T)\ 2 ~\u 2 (T)| 2 }, 

where HC stands for Hermitian conjugate, and 

ui(t) = —{ (e~ rt + e rt ) e* sin 9 + re"* 7 (e~ rt - e rt ) cos 9 } , 



u 2 (*) = £_{ ( e - rt + e rt ) cos 9 + re 1 ^ (e~ rt - e rt ) e* sin 9} , 
re* 7 = fcu + iA 12 . 

In this calculation, we used the parameter differentiation technic, explained in |f|. At time 
t = ti s + ^os + t 2s + T, atom A p , initially in state g, starts to interact with mode M\. 
Considering that A p interacts with M\ during a time t\ p = 3ir/ (2G), with M 2 during 
t 2p = {29 — 7r) / (2G), and a time t 0p = t 0s between these interactions, the probability to find 
the atom A p in state e at t — t\ s + t 0s + t 2s + T + t lp + t 0p + t 2p will be 

— 2kT 

P e = I (e~ rT + e rT ) + sin 29 cos (7 + <p) (e~ rT - e rT ) \ 2 . 

This is the main result of the present contribution. Observe that P e — 1 if there is no 
environment, since in this case k — r — 0. 



The probability P e depends on 9 and 4>, which may be freely chosen by varying the times 
tig, to s and t 2s . For fixed T and 6, a P e x <fi plot constructed with experimental data tells us 
about the cross decay rates: if any sinusoidal pattern is observed, k 12 + iA 12 is non zero, the 
amplitude of the curve being related to r; the position of the maximum of the curve may 
be used to find 7 (this maximum occurs for 7 + = tt). This is exemplified in Fig. |21 The 
larger visibility is achieved for 6 = 7r/4, i. e., when the state (J^J) is maximally entangled. 
Since the magnitudes of the cross decay rates are related to correlations between the way 
each mode interacts with the environment, we expect that the amplitudes of the curves grow 
when the cavities get closer. 

In order to get better insight into how the robust states arise in the model, let us rewrite 
Eq. © in the form 



in 
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C = d+C 2 , 

Ci = (k - r) (2Ai • A\ - •A[A 1 - AjAj 
C 2 = (k + r) (2A 2 • A| - *A\A 2 - A\A 2 

where the operators Ai and A 2 , given by 

II 1 -exp(-ry) 
exp (47) 1 

obey the usual commutation relations for bosons. Notice that the decay rate related to C± 
is smaller than the one related to C 2 . Thus, states which may be written as 
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(7) 



are less exposed to the environment. When r = k, these states are decoherence-free, a 
probably very difficult condition to reach. If we learn the value of 7, by performing the 
experiment just described, we may chose to work with the states less exposed to the en- 
vironment (jZJ). Relevant examples for Quantum Information and Quantum Optics are the 
maximally entangled state 



and the coherent state 



p s = (|l,0>-e^|0,l» ((1,01-e-^ (0,1|) 



Ps = \v, — e n v ) (v, — e ll v | 



As discussed above, the physical distance between the cavities may lead to the destruction 
of the microscopic correlations important for the appearing of DFS and robust states. There 
is, however, another possibility, which circumvents this problem: one can use two modes 
differing by their polarization in a single cavity. Is this a better proposition than the first 
one? In order to investigate this, let us consider two degenerate modes with orthogonal 
polarizations and equal dissipation rates. A circular Rydberg atom with e «-> g transition 
resonant with the modes will interact simultaneously with them via the Hamiltonian 

Hi2a = Kl (afai + a|a 2 + yj + ^ («a{er_ - iaia + + ai>cx_ + a 2 cr + j . 

The phase difference in the coupling constants of the atom to each mode are due to the 
polarization orthogonality. Consider that one atom in state e begins to interact with the 
modes, initially in vacuum state, at time t — 0, and this interaction lasts until the time 
t\2a = tt/(2\/2G). This atom creates the entangled state 

Ps (t = h 2a ) = (|0, 1} + % |1, 0)) (HC) (8) 

for the modes in the cavity, and ends up in state g. Taking into account the action of the 
environment in the period between t = tu a and t = ti 2(l + T ( T 3> n / (2\/2G)), the state 
of the fields will be given by Eq. (JBJ) with u\ (t) and w 2 (£) calculated using 9 = 7r/4 and 
cf) = tt/2. Then, if another atom, initially in state g, starts its interaction with the modes at 
time t = t\2a + T , and this interaction lasts until t = 2ti 2a + T, we get 

— 2kT 

Pe,r = 6 — | {e~ rT + e rT ) - sin ( 7 ) {e~ rT - e rT )\ 2 

for the probability to find this second atom in state e. Notice that P e = 1 if we disregard 
the environment (k = r = 0). Another way to have P e = 1 is k = r and 7 = 7r/2: in this 
case we have a DFS. 

If the cavity is subjected to a pressure, in such a way that it becomes slightly elliptical, 
the orthogonal modes become non resonant. For a large enough detuning, we may work 
with the Hamiltonian 



Hi2a = (a 2 a 2 + yj + HG (a^ 17 - + a 20"+) 



related to an atom interacting with just one resonant mode. Now, an atom, initially in state 
e, interacts with the field mode, initially in vacuum state, between t — and t = n/2G, 



yielding the field state |1) and ending up in state g. If we take into account the action 
of the environment between t = ir/2G and t = tt/2G + T (T ^> n/2G), and let a second 
atom, initially in g state, to interact with the field in the period between t = n/2G + T and 
t = ir/G + T, we get 

p - P -UT 
1 e,nr ° 

for the probability to find this atom in state e. 

The difference between probabilities D = P e>r — P e>nr must be a sign of significative 
cross decay rates. In the experiment with two cavities, we can maximize the effects of the 
parameter r by choosing the state built by the first atom; in this experiment with a single 
cavity, the first atom always constructs the state (JSJ), and the effects of r will be maximized 
for 7 = 7r/2 only. This is the condition that makes D more easy to detect, and it was used 
in Fig. El Since 7 depends on the environment, it is not trivial to control it. 

Cross decay is related to interference: if the environment acts over both modes in a 
microscopic correlated way, this action may be canceled (at least partially) for a set of states. 
For the ideal case of perfect correlation, the cross decay rates are large, and a DFS appears. 
If we get some correlation (even non perfect), the knowledge about the cross decay rates 
may be useful, since some states will be more resistent to the degradation produced by the 
environment. The experiments we proposed here are feasible with the present technological 
state, and the confirmation of these interference effects in Cavity Quantum Electrodynamics 
systems could encourage a search of analogous behavior in other contexts, maybe in scalable 
systems, which would be important for Quantum Information implementation. 
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FIG. 3: Difference between probabilities D = P eiT — P e ,nr for k = 1000 s" 1 , 7 = tt/2 and various 
values for r: r = 500 s -1 , r = 900 s" 1 and r = 1000 s _1 . Higher r values correspond to higher 
curves. 
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